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ABSTRACT 

Based on the study of temporal asymmetry of 631 gamma ray bursts from the BATSE 
3B catalog by Link and Epstein [Ap J 466, 764 (1996)], we identify the population of 
bursts whose rising times are longer than their decays, thus showing atypical profiles. 
We analyse their sky distribution, morphology, time-space clustering and other average 
properties and compare them with those associated with the bulk of the bursts. We 
show how most of the peculiar bursts analysed are consistent with recent fireball 
models, but a fraction of bursts (~ 4% of the total sample) appear to be inconsistent. 
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1 INTRODUCTION 

During more than 25 years, the origin of gamma ray bursts 
(GRBs) has been, perhaps, the deepest and most persistent 
problem in astrophysics. However, with the advent of the 
Compton Gamma Ray Observatory (CGRO) and its Burst 
and Transient Source Experiment (BATSE) in 1991, a new 
phase in the research of GRBs started. In seven years of 
operation, BATSE has accumulated a database of more than 
2000 observations. The angular distribution of these bursts 
is isotropic within the statistical limits, and the paucity of 
faint bursts implies that we are seeing to near the edge of 
the source population (e.g. Meegan et al. 1992, Fishman & 
Meegan 1995). Both effects, isotropy and non-homogeneity 
in the distribution, strongly suggest a cosmological origin of 
the phenomenon. In support of this conclusion, absorption 
lines (Fe II and Mg II) in the optical counterpart of GRB 
970508 have been detected with a redshift of z — 0.835. 
Along with the absence of Lyman-a forest features in the 
spectra, these results imply that the burst source is located 
at 0.835 < z < 2.3 (Metzger et al. 1997). 

The energy required to generate cosmological bursts is 
as high as 10 51 erg s _1 . The very short timescale observed 
in the time profiles indicate an extreme compactness that 
implies a source initially opaque (because of 77 pair cre- 
ation) to 7-rays. The radiation pressure on the optically 
thick source drives relativistic expansion, converting inter- 
nal energy into kinetic energy of the inflating shell. Baryonic 
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pollution in this expanding flow can trap the radiation until 
most of the initial energy has gone into bulk motion with 
Lorentz factors of F > 10 2 — 10 3 . The kinetic energy, how- 
ever, can be partially converted into heat when the shell 
collides with the interstellar medium or when shocks within 
the expanding source collide with one another. The random- 
ized energy can be then radiated by synchrotron radiation 
and inverse Compton scattering yielding non-thermal bursts 
with timescales of seconds. This fireball scenario has been 
developed by Cavallo and Rees (1978), Paczyhski (1986), 
Goodman (1986), Meszaros and Rees (1993), Meszaros, La- 
guna and Rees (1993) and others. A comprehensive review 
is presented by Meszaros (1997). 

The fireball model is a robust astrophysical scenario in- 
dependent of the mechanism assumed for the original energy 
release. A popular mechanism is the merger of two collapsed 
stars in a binary system, for instance, two neutron stars or 
a neutron star and a black hole (see Narayan, Paczynski & 
Piran 1992, and references threin), although other processes 
have been suggested (e.g. Usov 1992, Carter 1992, Melia and 
Faterzzo 1992, Woosley 1993). 

One important prediction of the fireball model, as well 
as by any explosive mechanism, is that individual burst pro- 
files should be inherently asymmetric under time reversal, 
with a shorter rise time than the subsequent decay time. 
This is a natural consequence of a sudden particle energy 
increase (e.g. produced by a shock) and the slower radiative 
dissipation of the energy excess. 

Time asymmetry in GRBs light curves has been dis- 
cussed by several authors (e.g. Mitrofanov et al. 1994, Link 
et. al. 1993, Nemiroff et al. 1994). In particular, Nemiroff et 
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al. (1994) showed that in the sample formed by those bursts 
with count rates greater than 1800 counts s _1 and durations 
longer than Is detected by BATSE until 1993 March 10, 
there is a significant asymmetry in the bursts profiles in the 
sense that most bursts rise in a shorter time than they decay, 
in agreement with what is expected from a general fireball 
model. The most recent and complete study was made by 
Link and Epstein (1996). They took 631 GRBs from the 
BATSE 3B catalog, including both faint and bright bursts, 
and confirmed the global asymmetry in the burst profiles 
showing that about two thirds of the events display fluxes 
that rise faster than the subsequent fall. About 30% of the 
bursts, however, presented a peculiar asymmetry in the tem- 
poral profiles, with slower risings than decays. 

In this paper we focus on this subsample of peculiar- 
asymmetric bursts (PABs) , which seems at first sight to con- 
flict with some predictions of the simplest scenarios for fire- 
balls. In particular, we shall discuss whether there are rea- 
sons to consider this subsample of GRBs as representative 
of a class of sources with different physical properties than 
other bursts. The structure of the paper is as follows. In Sec- 
tion 2 we define the sample and present the results of the 
symmetry analysis. We provide tables with the full results 
for PABs in order to allow identification of specific events. In 
Section 3 we study the sky distribution of the sample, while 
in Section 4 we investigate the level of positional coincidence 
(possible repetition) that PABs show. Finally, we discuss the 
implications of these results for theoretical models of GRBs. 



2 SAMPLE AND SYMMETRY ANALYSIS 

We have studied the sample of 631 bursts from the BATSE 
3B catalog whose global symmetry properties were discussed 
by Link and Epstein (1996). This sample contains both faint 
and bright bursts, spanning a 200-fold range in peak flux. 
PREB plus DISC data types at 64 ms time resolution, with 
four energy channels, were used in the analysis. 

The time asymmetry of the individual burst profiles was 
examined with the skewness function introduced by Link et 
al. (1993) and used in Link and Epstein's (1996) paper. This 
function is defined as 



A = 



<{t- <t >) 3 > 

<{t-<t >) 2 >3/2' 



(1) 



Here, angle brackets denote an average over the data sample, 
performed as 



< g(t) >5 



~ c th)g{u 



(2) 



where a is the measured number of counts in the ith bin, ti 
is the time of the ith bin and c t h is a threshold level defined 
as. 



cth = f{cp - b) + b. 



(3) 



Here c p stands for the peak (maximum) count rate, b is the 
background, and / < 1 is a fraction that will be fixed for 
the data set. Fixing / ensures that A is calculated to the 
same fraction of the peak flux relative to the background. 
Larger values of / emphasize the structure of the peak over 
the surrounding foothills. The normalization of A makes it 
independent of background, duration, and amplitude. It is 



equal to in the case of symmetric bursts, greater than 
for a burst whose peak rises more quickly than it falls and 
smaller than in the opposite case. It is equal to 2 for an 
exact FRED (from fast rise and exponential decay) and to 
—2 for an exact anti-FRED. Four fixed values of / were 
analyzed: /i=0.1, / 2 =0.2, / 3 =0.5 and / 4 =0.67. The only 
requirement for a burst to be tested in each of the four /- 
values is that the number of bins whose d exceeds or equals 
c t h is at least three. Consequently, the size of the sample 
differs for each choice of /. The error bars in A represent la 
deviations, calculated by randomizing the number of counts 
according to Poisson statistics and computing the variance 
of the asymmetry parameter for many trials. 

In Tables 1 - 4 we show the results of the skewness 
analysis for those bursts that presented PAB behavior (i.e. 
A < 0). Each table contains the peak flux, trigger number, 
burst type, value of the symmetry parameter for each fi, 
and a classification of the bursts profiles according to the 
following scheme: S for single-peaked or spike-like bursts, M 
for multiply-peaked bursts, and C for complex or chaotic 
events. Regarding the burst types, events for which the A is 
negative for all / are labelled "1" , whereas events for which 
the errors in A allow positive A for at least one value of / 
are denoted "2" . In Fig. 1 we show specific examples of these 
profiles. 

We found that 91 out of 631 bursts (14.4%) are PABs, 
i.e. do not present positive skewness for any /JT| Only 28.5% 
of these bursts are single-peaked. The rest are multiply- 
peaked or complex events. Notice that most of S-type bursts 
are in Table 4. This is consistent with the analysis tech- 
nique: a fast burst, typically lasting a couple of seconds, will 
have few points above the higher cut-offs and then, data for 
Af 2 j 3 ,f 4 will not be computed. 

As we can see from Tables 1 - 4 as well as from Fig. 
1, PABs exhibit a variety of temporal morphologies. If all 
of these events are produced by a single mechanism, then 
there should be a very wide range of boundary and initial 
conditions in the sources in order to generate such a plu- 
rality of profiles. With the aim of searching for differences 
between PABs and the more common bursts, we have com- 
puted average values of the hardness ratio and durations of 
type 1 PABs. These values are compared with similar esti- 
mates for those bursts with A > at all levels in Table 5. 
Due to the small number of bursts and the variety in their 
features, dispersions are so large that no conclusions can be 
drawn. However, in the light of current data, it is clear that 
no significant correlation is found between hardness ratio, 
or duration, with burst morphology. 



3 SKY DISTRIBUTION 

One of the most important results of BATSE is the discovery 
of that GRBs are isotropically distributed on the sky (see, 
however, Balazs et al. 1998). With the recent detection of 
high-redshifted absorption lines in the optical counterparts 
of individual bursts (Djorgovski et al. 1997, Metzger et al. 
1997, van Paradjis et al. 1997) Galactic models appear to be 
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finally ruled out. However, one could ask whether the distri- 
bution of PABs exhibits the same level of isotropy than that 
of the whole sample. It could be the case, for instance, that 
PABs have a different origin than other GRBs, and conse- 
quently, display a distinct distribution on the sky (e.g. there 
could be a statistically significant concentration of PABs in 
the supergalactic plane or within any superstructure). 

In order to quantify the isotropy we followed the method 
developed by Briggs (1993). The dipole moment toward the 
Galactic center is < cos 9 > , the mean of cos , where Qi is 
the angle between the ith burst and the Galactic center. An 
excessively large value of < cos# > indicates a significant 
dipole moment towards the Galactic center. The quantity 
< sin 2 b — 1/3 > tests for a concentration in the Galactic 
plane or in the Galactic poles. The expected mean values 
of the two statistics are zero for an isotropic distribution 
and, if they are asymptotically gaussian distributed, i.e. if 
for a large number of bursts in the sample (N) they are 
gaussian distributed, the variances a 2 are 1/3N and 4/45./V 
respectively. 

Briggs et al. (1996) noted that because the CGRO is in 
a low-Earth orbit, about one-third of the sky is blocked by 
the Earth causing a portion of the Galactic equator to be 
observed about 20% less than the poles. An additional effect 
is different exposure times between the Galactic south and 
north poles. These effects must be taken into account when 
computing the expected values of the statistics (Briggs et al. 
1996). Location errors on particular bursts, however, have 
no impact on the isotropy characteristics because they are 
small compared with the large scale of anisotropies we are 
testing against. 

Table 6 shows that the distribution of all PABs (Fig. 
2) is consistent with perfect isotropy. The same is true for 
sub-samples of PABs. Some entries in Table 6 show small 
deviations from isotropy (quadrupole). However, the small 
number of events make the asymptotic gaussian distribution 
no longer valid, and one should compare with the study of 
Briggs et al. (1996) (see their Fig. 4a and b). Comparing 
with the values of a that arise from the previous cited figures 
of the Briggs et al. work, we find, consequently, that there 
is no detectable anisotropy in the sky distribution of PABs 
and we see that the la deviation from isotropy contains the 
values of all entries in Table 6. 



4 TIME-SPACE CLUSTERING 

Several time-space clustering analysis of different GRB- 
samples have come to contradictory conclusions about 
whether some GRBs repeat or not (e.g. Quashnock & Lamb 
1993, Narayan & Piran 1993, Wang & Lingenfelter 1995, 
Petrosian & Efron 1995, Meegan et al. 1995). The most 
complete study on the subject until now, carried out by 
Tegmark et al. (1996), is based on the analysis of the angu- 
lar power spectrum of 1120 bursts from BATSE 3B catalog. 
These authors found that the number of bursts that can be 
labelled as repeaters (considering just one repetition) is not 
larger than 5% at 99% confidence. The recent study of by 
Gorosabel et al. (1998), which combined data from different 
satellites, shows that at most 15.8% of the events detected 
by WATCH recur in the BATSE sample (at 94% confidence 
level). Despite the discussion in the literature, it seems clear 



that only a small fraction of the total number of GRBs could 
repeat over timescales of up to a few years. 

However, if PABs have a different physical origin than 
other bursts, this subclass of bursts might exhibit time-space 
clustering. In fact, we find that 48 out of 91 PABs (52.7%) 
have companions within their location error boxes in the 
sample of 631 bursts. If we consider just bursts separated 
by less than 4°, we find 40 possible repeaters (44% of the 
PAB-subsample); typically, the separation is about 2.5°. 

To estimate the statistical significance of these results, 
we have made a numerical study as follows. We have sim- 
ulated 1500 sets of 91 random positions for PABs. In order 
to do this, we have made rotations on the celestial sphere 
sending a particular PAB with coordinates (I, b) to a new po- 
sition (l',b'), which is obtained from the previous by seting 
I' — I + Ri 360° and b' = b + R2 90°, and using appropriate 
spherical boundary conditions. Here, Ri and R2 are differ- 
ent random numbers (between and 1) which never repeat. 
Doing this for each event we get a new set of simulated 
PAB-positions. For this set we then compute the positional 
coincidence level with respect to the fixed 631 — 91 GRBs 
coordinates. As in the real case, we shall assign a positional 
coincidence when two or more bursts are separated by less 
than 4°. After making 1500 operations of this type (a larger 
number of simulations does not significantly modify the re- 
sults) we can obtain the mean value of the expected number 
of positional coincidences and its a. We obtain that for 91 
GRBs, the average level of positional coincidences is 42.9 ± 
4.7, which is entirely compatible with the observed result for 
PABs within ler. 

We have repeated the process for the subset of 26 single- 
peaked PABs (those denoted by an S in Tables 1 - 4) . These 
events represent about 4% of the whole sample and about 
28.5% of the PAB subsample. 15 out of 26 bursts of this 
kind (~60%) present companions within error boxes of less 
than 4 degrees. We find that the average simulated positional 
coincidence level is 13.3 ± 2.5. That is, the real coincidence 
level is also compatible with the random result to within lcr 
and no particular association appears obvious. 

If we now take positional coincidences separated by less 
than 1°, we find that 3 out of 26 single-peaked PABs have 
companions. Repeating the simulations in this case yields an 
expected chance association of 1 ± 1 events. This means that 
the real positional coincidence is only compatible with the 
random one to within 2a. The number of events is of course 
too scarce to draw any conclusion, but if this is confirmed in 
a larger sample it would entail an excess of 3.8% repetitions 
above the result expected from chance associations (some- 
thing compatible with already mentioned Tegmark et al.'s 
analyses) . 

As we shall see in the next section, spikes with peculiar 
asymmetry present problems for their interpretation within 
the standard fireballs models. 



5 PABS AND THE FIREBALL MODELS 

As mentioned in the introduction, GRBs profiles with A < 
are not expected from the simplest versions of the fireball 
model (i.e. a single expanding shell that acts as a gamma 
ray emitter during a brief time at some fixed radius from the 
central site of the explosion, e.g. Fenimore et al. 1996). How- 
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ever, one of the distinctive features of the fireball scenario 
is that the same basic mechanism can generate a variety of 
time profiles for different initial and boundary conditions. 
We now discuss whether these changes can provide the main 
types of PABs observed in the sample. 

In Fig. 3 we show the profile of BATSE trigger #2450, 
which has negative skewness function for all values of / (see 
Table 1). This is a typical multiply-peaked burst, with a 
precursor at t = and a series of peaks of increasing height 
that start about 35 s after the first signal. Individual peaks, 
when analysed with appropriate /, give A > 0. Events of this 
kind can be understood as the effect of a mild baryon loaded 
fireball (Meszaros & Rees 1993). Even a small baryon con- 
tamination (Mf, > W~ 9 Mq) of the expanding pair-photon 
fireball is enough to trap the 7-rays until most of the ini- 
tial energy is transformed into kinetic energy of the baryons. 
The fireball expands by radiation pressure and becomes op- 
tically thin to Thomson scattering when the optical depth 
drops below unity at a radius given by (Meszaros & Rees 
1993), 



0.6 x 10 1 



(4) 



where takes into account the possibility of channeling of 
the flow (9 ~ 1 corresponds to spherical symmetry), E51 is 
the original energy release (e , 7) in units of 10 51 erg, and 
r\ — Eo/Moc 2 is the initial radiation to the rest mass energy 
ratio in the fireball. At r — r p , the 7-rays still trapped in the 
fireball can escape producing a burst (Cavallo & Rees 1978, 
Paczhski 1986, Goodman 1986). As shown by Meszaros & 
Rees (1993), this burst should be rather weak, with an ob- 
served energy in gamma rays of, approximately, 



E° p bs ~ 7 x ia"0 1 '*E%{ a rg L/e erg, (5) 

where 773 = 10 -3 ??. This prompt, small burst will form a 
precursor that can last a few seconds. When the expand- 
ing relativistic shell collides with the interstellar medium, a 
shock wave is formed and the gas in the post-shock region 
is heated up to thermal Lorentz factors of 7 ~ 77, reconvert- 
ing the kinetic energy of the shell into thermal energy of 
the particles. The thermal energy is radiated through syn- 
chrotron and inverse Compton processes at MeV to GeV en- 
ergies. A non-uniform ambient medium can naturally lead to 
a multiply-peaked burst (e.g. Fenimore et al. 1996). Events 
of this class will have A < 0, as in the case of trigger #2450, 
due to the effect of the prompt precursor. Hence, these A < 
events can be explained within the fireball model. 

In other cases, the precursor can remain undetected but 
a multiply-peaked PAB can arise from internal shocks in 
bursts with several shells with different Lorentz factors (e.g. 
Kobayashi et al. 1997, Daigne & Mochkovitch 1998). In the 
simulations carried out by Kobayashi et al. (1997), bursts 
with negative skewness can be produced through multiple 
shell collisions (e.g. see Fig. 2f of their work). 

Complex bursts, as the one shown in Fig. Id, could be 
the result of instabilities on the expanding shell surface once 
it shocks the interstellar medium. Hydromagnetic instabili- 
ties in the contact discontinuity can lead to local variations 
in the fields and the flow's Lorentz factor, yielding very rapid 
changes in the time profiles (e.g. Daigne & Mochkovitch 
1998). The resulting global morphology could resemble that 
seen in some bursts with negative skewness, such as #2240. 

Single-peaked bursts with A < 0, however, appear to be 



more difficult to explain with the fireball model. The main 
problem is that a single spike with slower rising than falling 
cannot be generated through dissipative shocks. In Fig. 4 
we show BATSE trigger #444 (see also Table 3 and Fig. 
la). We have attempted to fit this event with the multiple 
shell model developed by Kobayashi et al. (1997). The 7-ray 
emission is produced when a shock results from the colli- 
sion of two shells with different velocities. The randomized 
kinetic energy is then radiated through synchrotron and in- 
verse Compton processes. Notice that the better the fit for 
the rising profile, the worse the model describes the fall. 
This is a straightforward consequence of the fact that cool- 
ing times are longer than particle acceleration times at the 
shock. 

To better understand the meaning of the theoretical 
curves in Fig. 4 we recall the predicted luminosity in the 
case of a two shell interaction (Kobayashi et al. 1997), 

C I - (l + 2jict/R)- 2 , 0<t<St e /2^ 
£(t)oc<^ (l + (2-y 2 n t-St e )c/R)- 2 -(l + 2-ylct/R)-' 2 , (6) 
{ t > <5i e /27m 

where 7™ is the Lorentz factor of the merged shell (depend- 
ing on the Lorentz factor and mass of each colliding shell), 
5t e /2'j 2 n is the time at which the burst reach its maximum, 
and R is the radius at which the collision takes place. Ob- 
servational data of a given burst, its height and duration up 
to the maximum in the number of counts, allow a parame- 
terization of C(t) with 



B = 



27m c 

R ' 



(7) 



The shape of the pulse is asymmetric with a fast rise and a 
slower decline unlike a spike event with A < 0. Attempts to 
fit such a burst using eq. (6) are shown in Fig. 4. 

Spike-like bursts with A < are predicted, however, 
in some extrinsic models for GRBs. Torres et al. (1998a,b) 
have shown that microlensing effects produced upon the core 
of high redshifted AGNs by compact extragalactic objects 
which violate the weak energy condition at a macroscopic 
level would yield GRB-like lightcurves with spike-type pro- 
files and negative skewness function. A similar burst with 
.4 > should be observed from several months up to a few 
years later in the same position of the sky, provided the lens 
has an absolute mass of the order of IMq . If this interpreta- 
tion turns out to be correct, it could explain not just S-type 
PABs but also any apparent excess of positional coincidences 
among these bursts at a level compatible with current con- 
straints on repetition over the whole sample. It would appear 
that the small group of spike bursts with negative skewness 
deserve further study. 



6 CONCLUSIONS 

GRBs exhibit a very rich variety of temporal profiles. Most 
of them have highly variable structure over timescales sig- 
nificantly shorter than the overall duration of the event. The 
study of burst morphology by Link and Epstein (1996) shows 
that a significant fraction of bursts (~ 1/3) have time his- 
tories in which the flux rises more rapidly than it decays 
(PABs). Here we have argued that most PABs can be ac- 
commodated by fireball models. Isotropy and other aver- 
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age features, common to the bulk of observed bursts, are Woosley S. 1993, Ap J 405, 273 
shared by PABs. But there is, however, a subclass of PABs, 
those which consist of a single, prominent peak with nega- 
tive skewness, that appear to be inconsistent with the fire- 
ball mechanism. These events represent ~ 4 % of the total 
sample and certainly merit further research in order to clar- 
ify their nature. 
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Table 1. GRBs with negative skewness function for all four cut-off parameters. Type 1 indicates bursts whose errors in A are less than 
the modulus of A itself and type 2 indicates cases in which at least one of the computed errors in A exceeds the modulus of A. S, M, and 
C stand for single, multiply-peaked and complex temporal profiles. 
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1.64800 


2074 


i 


-9.33E-02 ± 4.34E-02 


-0.19 ± 9.48E-02 


-0.31 ± 0.11 


-0.25 zh 0.20 


M 


1.54500 


254 


2 


-3.68E-02 ± 7.61E-02 


-6.36E-02 ± 8.59E-02 


-0.15 ± 0.19 


-0.19 ± 0.36 


S 


1.44900 


1197 


2 


-0.27 ± 4.82E-02 


-0.29 ± 0.14 


-0.17 ± 17 


-0.25 ± 0.27 


c 


1.37400 


1413 


2 


-0.18 ± 0.20 


-0.29 ± 0.24 


-2.68E-02 ± 0.30 


-0.14 zh 0.27 


M 


1.35300 


2610 


1 


-0.31 ± 3.31E-02 


-0.19 ± 3.28E-02 


-0.87 ± 6.81E-02 


-0.39 ± 0.21 


C 


1.11400 


752 


2 


-0.19 ± 0.14 


-0.24 ± 0.13 


-0.61 ± 0.29 


-0.32 ± 0.34 


s 


1.11000 


2828 


2 


-7.11E-02 ± 0.29 


-9.19E-02 ± 9.92E-02 


-0.14 ± 0.16 


-0.15 ± 0.27 


c 


0.84500 


2448 


2 


-9.77E-02 ± 1.87E-02 


-0.11 ± 2.28E-02 


-7.29E-03 ± 7.45E-02 


-0.11 ± 8.27E-02 


c 


0.77600 


2495 


1 


-2.43E-02 ± 1.23E-02 


-4.93E-02 ± 1.22E-02 


-0.15 ± 0.14 


-1.20 ± 0.68 


c 


0.66100 


690 


2 


-0.63 ± 0.21 


-0.35 ± 0..32 


-0.19 ± 0.37 


-0.30 ± 0.46 


c 


0.59400 


2442 


2 


-2.36E-02 ± 3.01E-02 


-1.03E-02 ± 0.11 


-0.13 ± 0.17 


-0.22 ± 0.14 


c 



Table 2. GRBs with negative skewness for three thresholds /. For f^, these bursts did not fulfill the criteria for the 
assignment of an A value. Burst type and capital letters have the same meaning as in Table 1. 



Peak flux 


# 


Type 


A/ x 


A/ 2 


A/3 


Profile 


12.45800 


1440 


1 


-7.47E-02 ± 8.56E-03 


-0.15 ± 8.27E-03 


-0.35 ± 0.17 


M 


8.73700 


2799 


1 


-0.78 ± 6.23E-03 


-1.1 ± 0.23 


-0.29 ± 8.58E-02 


M 


2.37600 


2861 


2 


-2.27E-02 ± 7.87E-02 


-0.12 ± 0.12 


-7.99E-02 ± 0.27 


M 


2.35100 


2795 


1 


-0.13 ± 1.62E-02 


-0.14 ± 1.92E-02 


-0.14 ± 5.31E-02 


S 


2.00500 


1359 


2 


-0.11 ± 0.14 


-0.14 ± 8.14E-02 


-0.76 ± 0.27 


S 


1.93900 


1154 


2 


-6.24E-02 ± 6.27E-02 


-8.17E-02 ± 6.22E-02 


-0.44 ± 0.24 


S 


1.23000 


1924 


1 


-0.22 ± 0.12 


-0.16 ± 0.11 


-0.17 ± 0.13 


S 


1.22600 


1611 


2 


-0.73 ± 0.34 


-0.20 ± 0.29 


-0.56 ± 0.26 


c 


1.20400 


3012 


2 


-2.87E-02 ± 0.12 


-0.15 ± 0.16 


-9.32E-02 ± 0.22 


M 


1.20200 


3080 


2 


-0.90 ± 0.36 


-0.17 ± 0.20 


-0.57 ± 0.25 


C 


1.19100 


2040 


2 


-4.77E-02 ± 0.16 


-5.50E-02 ± 0.18 


-8.62E-02 ± 0.25 


C 


1.05000 


3017 


2 


-6.30E-02 ± 0.10 


-0.17 ± 0.17 


-0.10 ± 0.19 


C 


0.88000 


2204 


2 


-9.25E-02 ± 0.28 


-0.41 ± 0.31 


-0.29 ± 0.22 


C 


0.81300 


2240 


2 


-0.18 ± 3.84E-02 


-0.20 ± 5.62E-02 


-0.10 ± 0.12 


C 


0.74000 


114 


2 


-0.13 ± 0.12 


-0.11 ± 0.28 


-0.38 ± 0.45 


C 


0.67800 


2434 


1 


-0.31 ± 2.74E-02 


-0.37 ± 3.15E-02 


-1.0 ± 0.24 


s 


0.54600 


2727 


2 


-7.14E-03 ± 0.11 


-1.43E-02 ± 8.57E-02 


-0.33 ± 0.50 


c 


0.34200 


1435 


2 


-0.15 ± 0.25 


-0.48 ± 0.23 


-0.65 ± 0.26 


c 
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Gamma Ray Bursts with peculiar temporal asymmetry 



Table 3. GRBs with A < for two thresholds. 



Peak flux 


# 


Type 


•4/1 


Af 2 


Profile 


28.55900 


444 


1 


-0.23 ± 1.66E-02 


-0.45 ± 3.51E-02 


s 


4 25500 


2993 


2 


-0 53 + 58 

'7 < "7 '7 1_ \7 . '7 17 


-9 09E-02 + 5 33E-02 

-7 < \j <y i 1 yj — i_ <7 • 77i - y ) 


M 


3.26900 


2041 


1 


-0.61 ± 8.86E-02 


-1.1 ± 0.36 


M 


2.22800 


2220 


2 


_n 44 + 46 

V ' < i i [_ V7 - i '7 


-1 66E-02 + 15 

1 . \J)\JX -. 11 _ 1_ 17 . 1 " 


S 


1 89300 


2201 


1 


-0 26 + 15 

'7 . Z. v ' 1_ \ ' . 1 • ' 


-0 32 + 7 65E-02 

y ) • - j z, i_ i . v ' • ' i - y t 


g 


1 85500 

1 . 1. ■ '7 • ' V ' V7 


2788 


2 


-4 68R-02 + 010 

i . '7 17 1 -. V7 — 1_ V7 . 1 V7 


-5.79E-02 ± 1.80E-02 


s 


1.22900 


1142 


1 


-0 83 -1- 39 

UtUU [_ V7 . '7 •/ 


-1.1 ± 0.42 


S 


1.20400 


1968 


1 


-0 26 + 19 


-0 29 + 14 


s 


1.07700 


171 


2 


-7 1 8E-02 + 26 


-9 49E-02 + 27 

CJ . L ± E7 J_J YJ £i _1_ YJ . Zj 1 


c 


1 06800 


2529 


1 


-0 40 + 5 1 5E-02 


-0 35 + 1 1 


M 


1 06000 

1 . \7 V ' V7 \ f \ f 


2853 


2 


-3 79E-02 + 9 67E-02 

*7 . 1 <J J_J YJ £1 _1_ i7 . V7 1 1 J YJ £j 


-8 85E-02 + 14 

L7 . * 1 '.7 1 -. 17 — 1_ YJ i _L rt 


c 


96700 

V 7 . v 7 V.' i YJYJ 


237 


2 


-0 11 +0 57 


-0 38 + 46 

17 . ■ J 7' [_ 17 . i '7 


c 


0.97900 


2800 


1 


-0.21 ± 0.12 


-0.27 ± 0.17 


c 


0.82500 


2290 


1 


-0.47 ± 4.08E-02 


-0.52 ± 4.27E-02 


c 


0.80800 


2129 


2 


-2.30E-02 ± 0.20 


-0.16 ± 0.23 


c 


0.78200 


1192 


2 


-2.31E-02 ± 0.13 


-3.93E-02 ± 0.20 


c 


0.78800 


2857 


1 


-0.25 ± 0.16 


-0.32 ± 0.18 


c 


0.76600 


1430 


2 


-8.19E-02 ± 0.26 


-0.10 ± 0.24 


c 


71600 

YJ • i 1UUU 


1404 


2 


-1.62E-02 ± 0.17 


-1.85E-02 ± 0.22 


c 


71 300 

YJ * t 1UUW 


1110 


2 


-0.11 ± 0.14 


-1.73E-02 ± 0.14 


c 


70800 

V 7 . 1 V 7 1. J 17 W 


204 


2 


-0 25 + 10 

UiiU [_ 17 . 1 V7 


-0.11 ± 0.14 


M 


70400 


2776 


I 


_n 23 + 8 03E-02 

U . Zj O _1_ O.UOJJ Uii 


_n 22 + 1 

U.ZjZi _1_ YJ . ±yj 


c 


65400 

17 . . 17 7 i 17 yj 


2996 


2 


-017 + 23 


-4 76E-02 + 48 

i . I vl 1 -. 17 — . 1_ V7 . i '7 


c 


65100 

17 . '7 7 1 '7 V7 


2725 


1 


-0 25 + 8 81-02 

'7 . £j "J 1_ O.tJJ-L V ' L- 


-0 40 + 19 

17 . i 17 i_ v7 . i v / 


c 


61800 

V 7 . 17 1 I. J V t\ / 


1655 


2 


-7 66E-02 + 31 

i .17111 j \ > ±- i_ v7 > v > i 


-0 37 + 28 

V 7 . - J 1 [_ 17 . — 17 


c 


0.57600 


2750 


2 


-4.14E-02 ± 0.13 


-0.11 ± 5.83E-02 


c 


0.56300 


2230 


2 


-0.10 ± 0.20 


-0.14 ± 0.18 


c 


0.56100 


2069 


2 


-9.99E-02 ± 0.15 


-0.37 ± 0.20 


c 


0.55200 


1382 


2 


-0.14 ± 0.18 


-0.28 ± 0.16 


c 


0.54500 


2900 


2 


-0.33 ± 0.24 


-0.12 ± 0.30 


c 


0.53400 


1465 


2 


-8.62E-02 ± 0.23 


-2 ± 0.34 


c 


0.49800 


2233 


1 


-0.48 ± 0.38 


-0.51 ± 0.13 


c 


0.47500 


559 


1 


-0.49 ± 0.31 


-0.31 ± 0.24 


c 


0.44000 


1120 


2 


-1.27E-02 ± 0.19 


-7.57E-03 ± 0.22 


c 


Table 4. GRBs with A < for one threshold. 






Peak flux 


# 


Type 




Profile 




8.76100 


2978 


1 


-1.60 ± 0.95 


S 




5.57300 


551 


1 


-9.78E-02 ± 5.33E-02 


s 




4.68100 


1851 


1 


-1.1 ± 0.44 


s 




4.08600 


2995 


1 


-0.83 ± 0.25 


s 




3.43700 


2918 


1 


-0.24 ± 0.20 


s 




2.97900 


297 


1 


-0.11 ± 3.65E-02 


M 




1.92600 


2161 


1 


-1.0 ± 0.63 


S 




1.79100 


2846 


1 


-0.36 ± 0.18 


S 




1.67000 


1461 


1 


-0.95 ± 0.41 


s 




1.13300 


2163 


2 


-0.13 ± 0.13 


s 




1.03200 


2823 


2 


-0.32 ± 0.34 


s 




0.69500 


2508 


1 


-1.6 ± 0.74 


c 




0.64800 


2437 


1 


-0.17 ± 4.50E-02 


c 





Table 5. Mean values of hardness ratio and temporal durations for type 1 PABs compared with usual A > bursts. All 
bursts have Af i computed for all four /; analysed. The PAB group has ~26 bursts, whereas the A > set comprises ~100 
of the total sample. 

<Hard. Ratio> a < T 90 > a 

type 1= PABs 3.6 2.3 35 42 

type 1= GRBs 3.2 1.5 49 50 
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Romero, Torres, Andruchow, Anchordoqui and Link 



Table 6. Isotropy characteristics. The expected values of the statistics (corrected for BATSE exposure) are -0.013 for 
< cos > and -0.005 for < sin 2 b - 1/3 >. 



Sample 


< cos > 


a = 0. 99^/1/3^3 


< sin 2 b - 1/3 > 


o- = 0.99A/4/457VS 


l+2-types=91 GRBs 


-0.034 


0.059 


-0.023 


0.030 


l-type=41 GRBs 


0.003 


0.089 


-0.045 


0.046 


2-type=49 GRBs 


-0.065 


0.082 


-0.005 


0.042 


1-type Table 1=12 GRBs 


-0.120 


0.165 


-0.052 


0.085 


1-type Table 3=13 GRBs 


0.072 


0.158 


-0.092 


0.081 


1-type Table 1+ 2=17 GRBs 


0.030 


0.138 


-0.104 


0.071 


1-type Table 1 + 2 + 3=30 GRBs 


-0.014 


0.104 


-0.098 


0.053 


single-peaked =26 GRBs 


0.015 


0.110 


0.015 


0.057 


Multiply-peaked and complex =65 GRBs 


-0.045 


0.071 


-0.044 


0.036 
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Gamma Ray Bursts with peculiar temporal asymmetry 



Figure 1. Examples of bursts with negative skewness. From top 
to bottom and left to right, Fig. la and lb are single-peaked 
bursts. Fig. lc is an example of a multiply-peaked burst with 
A < and Fig. Id is a complex peculiar asymmetric burst. For 
details on the particular yl-values see Tables 1-4. 



Figure 2. Spatial position (Galactic coordinates, Aitoff projec- 
tion) of the PABs presented in Tables 1-4. The sky distribution 
is consistent with perfect isotropy. 



Figure 3. An example of a burst with negative overall skew- 
ness, but positive skewness for the subpeaks. A burst of this type 
could be made up of a precursor associated with baryonic pollu- 
tion, followed by a series of peaks generated in the collision of a 
relativistic shell with the interstellar medium. 



Figure 4. An illustration of the inconsistency between the fireball 
model and spike bursts with negative skewness. Shown is the spike 
burst of Fig. la (see also Table 3). The bar plot stands for the 
observed number of counts. We also show several profiles of the 
light curves predicted by the multiple shell model of Kobayashi 
et al. (1997) corresponding to B = 0.1,5 and 15 (see text). In 
general, the better the theoretical curves fit the rising curve, the 
worse they fit the falling portion. 
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